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In astrophysical environments, the long-lived (T1/2 = 37.6 Gy) ground state of 176 Lu can com- 
municate with a short-lived (7i/ 2 = 3.664 h) isomeric level through thermal excitations. Thus, the 
lifetime of 176 Lu in an astrophysical environment can be quite different than in the laboratory. We 
examine the possibility that the rate of equilibration can be enhanced via jf-mixing of two levels 
near E x = 725 keV and estimate the relevant 7-decay rates. We use this result to illustrate the effect 
of Jf-mixing on the effective stellar half-life. We also present a network calculation that includes 
the equilibrating transitions allowed by Jf-mixing. Even a small amount of 7\-mixing will ensure 
that 176 Lu reaches at least a quasi-equilibrium during an s-process triggered by the 22 Ne neutron 
source. 

PACS numbers: PACS: 25.40.Lw, 26.20.-f 



I. INTRODUCTION 

The nucleus 176 Lu is produced almost exclusively by 
slow neutron capture in the s-process. It was originally 
thought that by virtue of its long half-life (T\/ 2 = 37.6 
Gy), 176 Lu could serve as a chronometer of galactic nu- 
cleosynthesis. However, it has been shown [j], Q that 
the ground state of 176 Lu can reach at least a quasi- 
equilibrium with a short-lived isomeric state at E x = 
122.9 keV (Ti/ 2 = 3.664 h; all excitation energies and 
other nuclear structure parameters are taken from the 
ENDSF compilation (3j) through a mediating level at 
E x = 838.6 keV. As a result, the astrophysical half-life of 
176 Lu is strongly dependent on temperature. This pre- 
cludes the use of 176 Lu as a chronometer, but does imply 
that it may serve instead as an s-process thermometer. 

The s-process in the in A — 176 mass region is shown 
in Figure [T] Production of the isomer, either by neutron 
capture or by thermal excitation of the ground state leads 
directly (via /3-decay) to 176 Hf. In contrast, the ground 
state, which is also produced directly by neutron cap- 
ture or by thermal excitation of the isomer, will undergo 
neutron capture to 177 Lu, which then /3-decays to 177 Hf . 
This temperature-dependent branching in the reaction 
flow should be reflected in the abundances in this mass 
region. In fact, consideration of isotopic abundances 
strongly suggests that the isomer and ground state must 
communicate during the s-process: As discussed by Klay 
et al. |1| (and references therein), measurements of the 
partial neutron capture cross section producing the iso- 
mer show that it is too large to be consistent with the 
observed abundance of 176 Lu in the ground state. Re- 
cent measurements reinforce this conclusion Stud- 
ies of the branching of the s-process at A = 176 also 



•present address: The University of Illinois, Urbana, IL 61801- 
3080, USA 



indicate that there is more s-process flow through the 
ground state than predicted by the production rate of 
the ground state @ via neutron capture in the classi- 
cal picture. Taken together, these observations suggest 
the presence of an additional creation mechanism for the 
ground state. Neutron capture produces an isomer to 
ground state ratio of approximately 5.7 [f|. At typical 
s-process temperatures (about 1-3 x 10 8 K, or Tg = 
0.1-0.3), the thermal equilibrium ratio is 1.3 x 10~ 7 -1.7 
x 10~ 3 . In other words, if thermal equilibrium pertains, 
then most of the s-process flow would ultimately pass 
through the ground state. 
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FIG. 1: The s-process path in the vicinity of 176 Lu. Stable (or 
extremely Hong-lived) nuclei are shown in grey. The ground 
state of 176 Lu and the isomeric state are labeled as g and m, 
respectively. 

Measurements of 7-decays by Klay et al. 0, Lesko 
et al. 2] and Vanhorenbeeck et al. Q have uncovered 
several paths linking the ground state and the isomer, 
which would be active during an s-process triggered by 
the 22 Ne neutron source. In this paper we will examine 
the possibility that iv-mixing provides a more efficient 
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way for the ground state and isomer to reach equilibrium. 
This mechanism is not limited to the case of 176 Lu and 
should operate in other s-process nuclei. 



II. EQUILIBRATION OF 176 Lu 

Direct 7-ray decay of the isomer (J v = 1 — , A 71 = _ ) 
to the ground state (J w = 7~, A 71 " = 7~) is strongly 
suppressed by their large difference in angular momen- 
tum. Also, because of the approximate conservation of 
the A quantum number, 7-decay is further hindered by 
a factor of 100 per degree of A-forbiddenness [§] (i.e., by 
a factor 10 2l/ , where v « \AK — A| and A is the 7-ray 
multipolarity). In general, the approach to equilibrium 
under stellar conditions can be inferred by comparing 
the time scales for 7-ray absorption and emission to the 
time scales for nuclear reactions and decays. States that 
are strongly linked by downward (and therefore upward) 
transitions, such as members of a given rotational band, 
will quickly reach full thermal equilibrium and will evolve 
together. States that are more weakly connected may not 
reach thermal equilibrium if other nuclear processes are 
faster. In the case of 176 Lu, the members of the ground- 
state band (tt7/2 [404]+^7/2 [514] , = 7") will quickly 
reach thermal equilibrium. For example, at Tg = 0.3, 
the 184.1-keV, = 8~ state will equilibrate with the 
ground state with a time constant of 21 ns. Side bands 
that feed the ground-state band will also equilibrate with 
the ground-state band. Similarly, the states in the isomer 
band (tt7/2 [404] - v7/2 [514] , K* = 0") and associated 
side bands will rapidly form an equilibrium cluster. At 
this point 176 Lu can be thought of as a quasi-equilibrium 
ensemble consisting of two separate, equilibrated clus- 
ters (the ground-state and isomer clusters). New groups 
of states in local equilibrium will form, presumably via 
A-allowed transitions (either upward or downward) from 
states in the original equilibrium clusters to those in other 
bands. If the time scale for linking transitions is short 
(compared to the time scale for neutron capture or (3- 
decay), then these new states will reach equilibrium with 
the original cluster. Otherwise, they will form separate 
groups in local equilibrium. Eventually, the clusters orig- 
inating from the ground state will reach those based on 
the isomer and produce quasi-equilibrium between the 
ground state and isomer. This process will be discussed 
in more detail in Sec. V. 

The previous 7-decay studies of 176 Lu, 0, 0] clearly 
show that the J* = 5~ state at E x = 838.6 keV (K = 4) 
decays directly to the ground state via a strong E2 tran- 
sition and via cascades that ultimately end in the isomer 
(more recent measurements [1[ show evidence for medi- 
ating states at somewhat higher excitation energy). In 
a stellar plasma, this state would be populated from the 
isomer in the manner described above, via 3-4 interband 
transitions, and primarily via a single photon from the 
ground state. Interestingly, several of the strongest tran- 
sitions from this level (including the ground-state branch) 



are A"- forbidden [l(|, which implies that the 838.6-keV 
state may be mixed with other nearby 5~ states. 

Klay et al. [l[ discuss the possibility that weak, un- 
detected A"-forbidden transitions may also provide an 
equilibration path, but argue that these contributions 
are small compared to the observed route. However, the 
effect of any such transition could be enhanced by Ad- 
mixing, which would largely circumvent the A-selection 
rule. As mentioned above, this mechanism appears to 
be at work for equilibration through the 838.6-keV state 
and there is clear evidence for A-mixing at higher excita- 
tion energies [ill ] . Although the lower-energy members of 
the isomer band have ratios of dynamic to kinematic mo- 
ments of inertia that are nearly unity (indicating little ro- 
tational perturbation), chance degeneracies can also lead 
to A-mixing and thus it seems at least plausible that the 
7- states at E x = 724.7 keV (tt7/2 [404]-i/7/2 [514] ,K* = 
Q~) and 725.2 keV (tt5/2 [402] -1/7/2 [514] , A T = 6~) will 
mix. The former state is a member of the isomer band 
whereas the latter state is a member of a band that can 
be populated by thermal excitation of the ground state. 
If these states do indeed mix, then the ground state and 
isomer are linked by a series of A-allowed transitions. 
As pointed out by Kondev et al. [H| and McGoram et 
al. [Ill ], even small amounts of mixing can significantly 
increase the rate of transitions that would otherwise be 
A-hindered. These transitions could be much weaker 
than the unmixed transitions and thus would have es- 
caped notice in previous decay studies. However, since 
the mediating states would be at lower excitation en- 
ergy, «725 keV rather than at 839 keV, they might be 
produced more efficiently from the stellar photon bath. 
This scenario for equilibration is illustrated in Figure [U 

III. A-MIXING OF 2 LEVELS AND 7-DECAY 

If we assume that 2 states of a given J 7r , but different A 
(Ai and A2) are mixed, then the resultant wavefunction 
for each state can be written schematically as 

\j*,Ki) =a\r,K 1 )-p\r,K 2 ) (1) 

and 

\j*,K 2 )=a\J",K 1 )+p\J",K 2 }, (2) 

where 

a 2 + /3 2 = l. (3) 

In our case, — 7~ ; Ai and K 2 refer to the mixed- 
A nature of the states (where A is no longer a good 
quantum number); the pure A states have K\ = 6, and 
A 2 = 0. 

The matrix element for the E2 transition linking e.g., 
states (5) and (6) in Figure[5] [henceforth, transition a— *b 
will be denoted as (ab)] is 

(5 - ,0||Af(£2)|| 7~,Ki) =a(5 - ,0||M(J52)|| 7~, 6) 

-/?(5-,0||Af(A2)||7-,0) (4) 
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FIG. 2: Proposed equilibration scheme. States are listed according to their J 71 " and excitation energies (in keV). The numbers 
in parentheses denote the states of interest. Downward 7-transitions have been observed for the links represented by black 
arrows. The matrix element that describes the mixing of the 7~ states [labeled (5) and (7)] in the K n = 6— and = 0— 
bands is denoted by V . The 2 bands are potentially linked by 7-transitions between states (5) and (6) and between (3) and (7) 
(shown in grey). 



where M(E2) is the electric quadrupole operator. The 
first term represents a i^-forbidden decay, which can be 
ignored. Assuming that T^ /2 /T^ 2 = lOfT (where T^ 2 
is the Weisskopf estimate for the partial half-life and v 
is the K- hindrance factor introduced above), T^ 2 (56) = 
488 ms, which is much longer than for the allowed decays 
calculated below. Since 



B{E2; J, K -> J -2,K) 
1 



\{J -2,K\\M(E2)\\ J,K)\\ 



(5) 



2J+ 1 
we have 

B(E2;7~,Ki -► 5 - ,0) = 1 B(E2; 7~ , -> 5",0). (6) 

From the relationship between B(E2) and the partial 
half-life for the transition, we finally obtain 



27/2(56) 



£ 7 (76)\ 5 I? /a (76) T7 /2 (76) 



£; 7 (56)y (P/a) 2 {(3 /ay 



(7) 



In other words, the i^-mixed (56) transition can be cal- 
culated from the partial half-life for the observed (76) 
transition if the mixing amplitude (3 can be determined. 

A similar approach can be used for the transition be- 
tween levels (7) and (3). However one difference is that 
the (73) and (53) transitions can be of mixed Ml + E2 



character. In this case 

1 1 



T i/ 2 ( 73 ) T i,2 ( M1 ; 73 ) T i,2 (£ 2 ; 73) 
_ l + <5 2 ( 73) 
" T7 /2 (M1;73)' 



(8) 



where 5(73) is the E2/M1 mixing ratio. Assuming that 
5(73) = 5(53), then 



37 /a (73) 



£ 7 (53)\ 5 I? /a (53) T7 /2 (53) 



E 7 (73)J (P/a) 2 {fi/af 



(9) 



IV. ESTIMATED 7-DECAY RATES 

The partial half-lives for the observed (53) and (76) 
transitions have not been measured and so we have calcu- 
lated them using a semi-empirical approach. For the (76) 
E2 transition, we use the relationship between B(E2) 
and the intrinsic quadrupole moment Qq: 

B(E2; J, K J - 2, K) = r|-Qg \(JK20\J - 2K)\ 2 

Ion 

(10) 

where Qo — 7.35(4) eb is calculated using the measured 
Q = 4.92(3) eb |3j for the ground-state band and the 
relationship 



3if 2 -J(J+l) 
W (J+l)(2J + 3)^°- 



(11) 
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The remaining term in Eqn. [TO] is a Clebsch-Gordan co- 
efficient. This yields B(E2) = 1.734 eb 2 , which implies 
T 7 /2 (76) = 16.7/a 2 ps. For the (53) Ml transition, we 
make use of 

B(M1; J,K -> J-1,K) = 

(g K -9R) 2 K 2 \(JK0\J~lK)\ 2 . (12) 

The gyromagnetic ratios §k (7r5/2[402]) = 1.57 and 
gx (^7/2[514]) = 0.30 were calculated using deformed 
Woods-Saxon wave functions with deformation param- 
eters /? 2 = 0.278, /? 4 = -0.071, and /3 6 = -0.009 The 
collective gyromagnetic ratio gn was calculated from 



A* = 9rJ + (9k - 9r) 



K 2 



J + l 



(13) 



with n = +3.1692(45) fj, N and \(g K - 9r)\ = 
0.175(7) [3. The resulting value, g R = 0.30 was as- 
sumed to be the same for all bands. From this, we obtain 
B(M1) = 0.298 n 2 N and T 7 /2 (M1) = 31.5 ps. However, 
since this transition can be of mixed Ml + E2 multipo- 
larity, 



<-pi 

J l/2 



T 7 

J l/2 



(Ml) 



1 + 5(53) 2 



where 



(5(53) = 0.933 



Qo £J 7 (53)(MeV) 



(9 



K 



VT 2 



1 



(14) 



= 0.3. (15) 



Therefore, T 7 /2 (53) 



9R) 

28.9/a 2 ps. 

The partial half-lives for the fC-mixed (56) and (73) 
transitions can be calculated using Eqns. [7] and [5] The 
mixing matrix element V and the amplitudes a and /3 
are related by 



V = a x /? x AE X , 



(16) 



where AE X is the energy difference between the two 7 
levels (517 eV). Given that a 2 + (3 2 = 1, Eqn.QUcan be 
solved for 1 : 



P 2 = 



2 ( 



1-4 



(A) 



1/2' 



(17) 



Unfortunately, V is unknown and can take values up to 
AE x /2 = 258.5 eV. Consequently, we have calculated 
partial half-lives and branching ratios for different values 
of V, which are listed in Table U 



V. EFFECTIVE DECAY RATE IN A STELLAR 
PLASMA 

The time evolution of the states shown in Figure [2] 
can be calculated by solving the set of coupled differen- 
tial equations that describe the behavior of each state. 



TABLE I: Mixing matrix elements and partial half-lives for 
A'-mixed transitions 



V 
(eV) 



P 2 
(%) 



T 



1/2 (76) 
(Ps) 



T 



1/2 (56) 
(Ps) 



T 



1/2(53) 
(ps) 



T 



1/2(73) 
(Ps) 



258.5 
150 
100 

50 

10 

5 

1 



50.0 

9.28 

3.89 

0.944 
3.74e-02 
9.40e-03 
3.74e-04 



3.34e+l 
1.84e+l 
1.74e+l 
1.69e+l 
1.671e+l 
1.6702e+l 
1.6700e+01 



3.34e+l 
1.80e+02 
4.29e+02 
1.77e+03 
4.47e+04 
1.79e+05 
4.47e+06 



5.78e+01 
3.19e+01 
3.01e+01 
2.92e+01 
2.89e+01 
2.8903e+01 
2.8900e+01 



5.78e+01 
3.11e+02 
7.43e+02 
3.06c+03 
7.73e+04 
3.09e+05 
7.73e+06 



Since various transition rates can differ by many orders 
of magnitude at s-process temperatures, significant care 
must be exercised to control the size of numerical errors. 
This becomes cumbersome once 176 Lu is placed into a 
larger network of nuclei. However, in a stellar plasma, 
176 Lu will rapidly evolve into a configuration of 2 quasi- 
equilibrium clusters, one of which contains the states 
in the ground-state and K — 6 bands, while the other 
is composed of the states in the isomer band. Within 
the clusters, the states rapidly achieve thermal equilib- 
rium, as described in Sec. II. To illustrate this point, we 
have calculated the time evolution of the states shown 
in Figure [5] for Tg = 0.25 and mixing amplitude 1 — 
9.40 x 10~ 3 %. Only the 184.1-keV and 388.9-keV states 
have tabulated lifetimes *3j and so those for the remain- 
ing states were calculated using the procedures described 
above. The results are shown in Figure [31 assuming that 
only the ground state is populated at t = and that 
there is no destruction of 176 Lu via /3-decay. The top 
panel shows the states that equilibrate quickly with the 
ground state. This cluster reaches thermal equilibrium 
within about 0.1 ns. The bottom panel shows the cluster 
built on the isomer, which reaches equilibrium somewhat 
later than the ground-state cluster (here the steady-state 
abundances are reached after about 1 ns) . The rise in 
the abundances of these states is a consequence of leak- 
age out of the ground-state cluster. At about 10 5 s, the 
isomer cluster reaches equilibrium with the ground-state 
cluster and all of the relative abundances assume their 
steady-state values. This time scale is about 14-15 or- 
ders of magnitude longer than the time scales for equi- 
libration within the clusters. Thus, on the time scale of 
equilibration between the ground state and the isomer, 
equilibration within a cluster can be considered to be in- 
stantaneous. Therefore, the only 7-ray transitions that 
have to be considered explicitly are those between the 
clusters, i.e. 5<->6 and 3<->7. 

In this approximation, the time evolution of the clus- 
ters can be described by 2 coupled equations: 




rii J = -ni\p - 713A37 - n 5 A 56 

+n 6 \65 + n 7 X 73 , (18) 
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FIG. 3: The time evolution of the states shown in Figure [2] 
for T 9 = 0.25 and (3 2 = 9.40 x 10~ 3 %, assuming no destruc- 
tion of 176 Lu via /3-decay.The top panel shows states in the 
ground-state cluster, labeled by their excitation energies (in 
keV). The isomer cluster is shown in the bottom panel. The 
two clusters equilibrate after about 10 5 s, which is 14-15 or- 
ders of magnitude longer than the equilibrium times within 
the clusters. The dots show the abundance of the isomer cal- 
culated using the 2-level approximation described in the text. 



which describes the ground-state cluster and 




-n 2 \ 2 p - uqX, 



n 7 A 73 
n 3 A 37 + n 5 A 56 



for the isomer cluster. Here rii refers to the abundance 
of state i in Figure [21 The subscripts g and m denote the 
equilibrium clusters built on the ground state and isomer, 
respectively; Ay are the i — ► j transition rates, and Ai 
and A| are the laboratory /3-decay rates for the ground 
state and isomer. Although the /3-decay rates depend on 
the degree of ionization, the temperatures and densities 
of the s-process are low enough that this effect may be 
neglected. The transition rate linking an initial state 
i with a higher-lying final state / via absorption of a 
photon is 



A 



if 



9f 



A 



fi 



g. e E tf /kT _ 1 ' 



(20) 



where g = 2 J + 1, is the rate for spontaneous decay 



from / to i, E^f is the energy difference between i and 
/, and k is Boltzmann's constant. The corresponding 
total decay rate between i and / (including stimulated 
emission) is 



A 



A 



if 



fi 



1 



-E, f /kT 



(21) 



Since the states within a cluster are in thermal equilib- 
rium, their relative abundances are given by 



9i - 
dj = — = — e 
n, g 3 



-E*j)/kT 



(22) 



Therefore, the abundances appearing in Eqns. [T8l and [T9l 
can be written as fractions of the ground state or isomer 
abundances. This allows these equations to be recast as 
a concise 2x2 coupled system of equations for the abun- 
dances of the ground state and isomer: 



dt 



dt 



rii = auTii + ai 2 n 2 



n 2 = a 2 \ni + a 22 n 2 



or equivalently as a matrix equation 



d 



— n = A ■ n. 

dt 



The transition matrix elements 



(23) 



(24) 



(25) 



an 



a 2 i 



(— ^^~C , 5iA 56 — C31 A37) 
(I+C31+C51 + ...) 
(C51 A56+C31 A37) 

1 + (C42 + C 6 2+C 7 2 + ...) 



a 12 



a-22 



(C72A73+C62A65) 
(1+C31+C51+...) 

( — A^ — C72A73 — C62A65 



1 + (C 4 2 + C 6 2+C72 + ...) ' 

(26) 

where a\ 2 and 021 represent the effective m — > g and 
g — > m rates, respectively. The denominators in these 
expressions are formally a sum over all of the Cy for 
each cluster. In practice, only the lowest-energy states 
contribute to this sum and for s-process temperatures, 



(19) Cij < 1. Eq.H5]has the solution 



n(t) = Nie 



-Ait 



N 2 e 



-A 2 t 



(27) 



where and -Aj are the eigenvectors (normalized to 
the correct initial abundances) and eigenvalues of A, re- 
spectively. The smallest eigenvalue corresponds to the 
effective decay rate of the (quasi) equilibrated nucleus. 
The other eigenvalue (which is approximately equal to 
022) is the inverse of the time constant for quasi or full 
equilibration. The resulting effective half-life of 176 Lu is 
shown in Figured! We have also calculated the effective 
half-life for equilibration through the 838.6-keV state, us- 
ing r = 10-300 ps for this state [HI, which is consistent 
with what is reported by Klay et al. [l[ • As shown in Fig- 
ure [31 this two-level approximation accurately accounts 
for the evolution of the isomer. 
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FIG. 4: The effective half-life of the ground state of Lu as a function of temperature. The results for A-mixing, assuming 
fractions of 3.74 x 10 -4 %, 9.40 x 10" 3 %, and 50%, are represented by the solid lines. The shaded area is the effective 
half-life if 176 Lu equilibrates via the 838.6-keV state. Full thermal equilibrium is denoted by the dashed line. 



VI. ASTROPHYSICAL ASPECTS 

Assuming that the if -mixing scenario is in fact valid, 
our calculation shows that the onset of quasi-equilibrium 
and the transition to full equilibrium may occur at lower 
temperatures than previously thought [1]. In particular, 
176 Lu will reach full thermal equilibrium for Tg > 0.22- 
0.30 for the range of mixing considered here. In contrast 
full equilibrium through the 838.6-keV state occurs for 
Tg w 0.29-0.32. This is not too surprising since we as- 
sume that the mediating level is at a lower excitation 
energy: 725 keV compared to 839 keV. However, it is in- 
teresting to note that a mixing amplitude of only 9.40 x 
10 _3 % is sufficient to compete with equilibration via the 
838.6-keV state. Since the corresponding mixing matrix 
element of 5 eV is on the order of atomic interactions, 
it is fair to wonder if any nuclear state can be consid- 
ered pure. In fact, the recent observation of AK = 13 
mixing in 1 Lu with V « 19 eV suggests a minimum V 
in the range of 10-20 eV [15j. This general conclusion 
likely applies to other if -isomers in the s-process (e.g. 
180 Ta 16]). In other words, the degree of if- mixing that 
exists in rotational nuclei is probably sufficient to allow 
if -isomers to reach at least quasi-equilibria with their 
ground states during some phases of the s-process. In 
the case of 176 Lu, quasi- or full equilibrium is reached for 
temperatures characteristic of the thermal-pulse phase in 
low-mass AGB stars [l?} or the weak s-process in mas- 
sive stars [18| . In both cases, the major neutron source 



is the 22 Ne(a,n) 25 Mg reaction. 

The effective half-life shown in Figure 0] gives an indi- 
cation of the temperature range where equilibration will 
occur, but the actual behavior of 176 Lu in a stellar en- 
vironment will depend on the competition between the 
time scale for equilibration, r eq (the inverse of the largest 
eigenvalue in Eqn. 127]) and those for neutron capture, r„ 
and /3-decay, rp. If r eq is smaller than both t„ and 173, 
then 176 Lu will maintain thermal equilibrium with an ef- 
fective half-life as shown in Figured]. Otherwise, 176 Lu 
will achieve a quasi-equilibrium in which the relative pop- 
ulation of excited states can be quite different from the 
equilibrium values. In this case, the approach to equilib- 
rium will lag behind abundance changes from n-capture 
and /3-decay and it is necessary to include the equili- 
brating transitions (ai2 and 021 from Eqns. [23] and I23p 
explicitly when calculating the time evolution of 176 Lu. 
During thermal pulses in low-mass AGB stars (with peak 
temperatures of Tg = 0.25-0.35 and neutron densities of 
8 x 10 6 - 3 x 10 10 cm -3 ), the question of whether full or 
quasi-equilibrium pertains depends strongly on the am- 
plitude for A-mixing. This is shown in Figure [5] For 
Tg= 0.35, 176 Lu will be in full thermal equilibrium for 
the entire range of A-mixing considered here because r eq 
<C t„ and rp (here we have assumed a neutron density 
of 3 x 10 10 cm -3 ). However, this is not the case in the 
limit of weak A-mixing and lower temperatures. At Tg 
= 0.25, r eq « T,3 for 1 = 9.40 x 10 -3 % and a quasi- 
equilibrium will ensue. A similar situation pertains for 



7 



equilibration through the 838.6-keV state. For T < 0.2, 
T eq > T /3 even in the case of maximal mixing and thus 
176 Lu can not be ascribed an effective decay rate. 
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FIG. 5: The time evolution of the isomer and ground state 
- assuming no production or destruction - for T9 = 0.25 and 
0.35 and ^-mixing amplitudes of 9.40 x 10 _3 % (dotted lines) 
and 50% (solid lines). The time scales r e? , r n and rp are 
denoted by vertical lines. For Tg = 0.35, r eq is much shorter 
than either r n or 773, which indicates that the equilibrating 
transitions are fast enough to maintain thermal equilibrium 
in the presence of n-capture and /3-decay. However, this is 
not the case for 9.40 x 10~ 3 % mixing at T 9 = 0.25. 

To show schematically how the abundance Lu will 
evolve during the s-process we have performed post- 
processing calculations for a 3 M Q AGB star of solar 
composition during the thermal-pulse phase. The tem- 
perature and density profiles at the base of the convec- 
tive shell as well as the abundances of the major isotopes 
were calculated using the Monash version of the Mount 
Stromlo stellar-structure code (Ref. [13] and references 
therein) . The initial abundances of the s-process isotopes 
were calculated for temperatures and densities appropri- 
ate to the interpulse phase, during which the 13 C(a,n) 16 
reaction is the primary source of neutrons. Neutron- 
capture rates in the Yb-Hf region were taken from the 
KADONIS compilation [f|. Since both 176 Hf and 176 Lu 
are produced almost entirely by the s-process (there is 
a small p-process contribution to 176 Hf, which we will 
ignore), the ratio 176 Hf/ 176 Lu is a useful s-process diag- 
nostic because its value at the time of the formation of 



the solar system is reasonably well established. We have 
adopted the value 7.67(18), which is obtained from the 
tabulated abundances of 176 Lu and 176 Hf [23], corrected 
for the radioactive decay of 176 Lu (with T x / 2 = 37.6 Gy) 
over the past 4.55 x 10 9 years. The results of our cal- 
culation for the 9 th thermal pulse are shown in Figure [6] 
Here, the peak temperature is Tg ~ 0.24 and thus the 
22 Ne neutron source is only marginally activated. The 
neutron burst lasts for about 12 years, which is about 
10% of the duration of the entire convective phase. These 
results are not meant to represent a quantitative depic- 
tion of the s-process for this star. We have calculated the 
s-process for conditions found only at the base of the con- 
vective shell and have not included mixing. Our purpose 
is simply to show how 176 Hf/ 176 Lu behaves for differ- 
ent degrees of equilibration (a more realistic treatment 
of nucleosynthesis can be found in Ref. [ll[). The effec- 
tive decay rate does a reasonable job of describing the 
evolution of 176 Lu for the case of maximal mixing (50%), 
but only after about 4 y beyond the start of the burst. 
For the (relatively) low temperatures in this pulse the use 
of an effective decay rate does not accurately describe the 
behavior of 176 Lu for lower mixing fractions and thus, the 
transitions linking the ground-state and isomer clusters 
must be used instead. One thing that is clear in Figure [6] 
is that 176 Hf/ 176 Lu decreases as mixing becomes more 
efficient. The same is true if the peak temperature is in- 
creased. The reason for this is that the equilibration time 
drops as the mixing fraction or temperature is increased, 
which favors production of the ground state of 1 Lu. 
High levels of if-mixing appear to be ruled out because 
they lead to too much 176 Lu relative to 176 Hf. As pointed 
out above, minimal amounts of if -mixing within the sce- 
nario outlined here are sufficient to compete with the 
observed transitions via the 838.6-keV state. The actual 
amount of mixing required depends on temperature, but 
a mixing fraction on the order of 3.74 x 10 _4 % - 0.015% 
or equivalently, a mixing matrix element of only 1-6.6 eV 
is sufficient throughout the range of temperatures for the 
22 Ne neutron source. As shown in Table HJ the branching 
ratios for the relevant 7-ray transitions are probably too 
small to be observed, which complicates the use of 176 Lu 
as an s-process thermometer. The small value required 
for the matrix clement implies that if -mixing should at 
least contribute to the equilibration of 176 Lu. However, 
without experimental verification, it may be difficult to 
predict the amount of 176 Lu produced, or 176 Hf/ 176 Lu 
for a given temperature regime. On the other hand, if 
small amounts of if-mixing can have a dramatic impact 
on equilibration, then it seems plausible that this mech- 
anism will operate in other places in the s-process. 
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